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INTRODUCTION 



In 1953, H. H. Bleich published a paper in the "Journal of Applied 
Mechanics" entitled "Non-Linear Distribution of Bending Stresses Due to 
Distortion of the Cross Section". In this paper he derived a viable 

analytical solution to the problem of hull-deckhouse interaction. 

0 

Basically, he considered the hull and the deckhouse as separate 
beams which are forced to act together by shearing forces and by verti- 
cal forces resisting relative displacements of the two beams. The case 
of constant cross section of the beam is treated, and it is assumed 
that Navier T s hypothesis is valid for the deckhouse and hull separately. 
For two types of loading he considered (constant moment loading and 
equally distributed load) , solutions were in a qualitative agreement 
with the test results at midship. As one moves away from amidships or 
the center of the deckhouse structure, solutions were departing from 
the reality. 

It is proposed that including shear effects into Bleich ! s original 
two-beam deckhouse theory, it is possible to improve it, particularly at 
deckhouse ends. In order to confirm or disaprove this hypothesis, new 
theory which includes shear effects is developed, and the results are 
compared at midship and at the end of deckhouse. 

The theorem of Minimum Potential Energy is used in the analysis. 

The variational procedure established two coupled six order differential 
equation systems and the natural boundary conditions. The boundary value 
problem is solved for the displacement and then the stress distributions. 
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FUNDAMENTAL ASSUMPTIONS 



The two components of the combination of hull and deckhouse are 
treated as if each were a Navier beam which are forced to act together 
by shearing forces and by vertical forces. 

It is assumed that the stiffness of bulkheads or deck beams 
resisting relative vertical displacement of* the deckhouse is constant 
for the full length of the deckhouse, the magnitude of the stiffness 
being given by a spring constant K. In reality, the deck stiffness K 
vary along the deckhouse length due to the presence of structural 
bulkheads . 

The possibility of having different materials for hull and deck- 
house is not considered in this analysis. 

Shear deformation is accounted for in the sides of both hull 
and deckhouse. For this purpose, the following equations are employed. 



M i - - E + 

lw 



M 2 = - E + G> 

2w 



( 1 ) 

( 2 ) 



The longitudinal stress in the deck at the junction and the 
longitudinal stress in the deck-edge may differ. Bleich’s theory 
neglects this. To include the shear-lag effect in the deck-edge and 
deckhouse side, the following assumption is used, 



a i r a 2 



( 3 ) 



where r is the ratio of longitudinal deck stress at the junction to 
the stress at deck-edge. 
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Precise determination of "r" would require two-dimensional 
elastic analysis of the response of all plate elements of the section. 

For the present purpose, it is assumed that "r" has the same value it 
would have without the deckhouse, determined by a box-girder analysis 
of the hull alone, such as in Reference (2). According to this analysis, 
for a sinusoidal bending moment: 

r = i ~ib ( Y sinh it + 2 cosh Y ■ ^ tanh cosh (4) 

cosh — 

It is pointed out in Reference (4) by Shade, for a bending 
moment which is constant over the length of the deckhouse "r" should 



be taken as unity. 
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CHAPTER I - ANALYSIS OF A SIMPLIFIED TWO CELL STRUCTURE 



Consider the problem of two separate beams forced to act together 

by horizontal shear forces and vertical forces acting at the junction of 

hull and deckhouse, (Figures 1 and 2). The vertical forces are due to 

elastic resistance of the deck framing or bulkheads against the motion 

* 

of the super structure with respect to the hull. The system consists, 
therefore, essentially of a beam elastically supported by another beam, 
with a shear connection to enforce equal strains at deck level. 

In this section the important simplifying assumption made is that 
the deck A-B, Figure 1, and its supports have no stiffness, and will not 
resist any relative vertical movements between hull and deckhouse. 

This simplification is not justified for any real ship system, but 
because of its relative simplicity it is easier to study the play of 
forces; the understanding gained is of value in treating the full 
problem in the following chapters. 

In the structure shown in Figure 1, the lower hollow box beam 
represents the hull and is of length L while the upper box, the 
deckhouse, is shorter and is of length £. Both boxes are assumed to 
be of constant cross-section. The cross-sectional area and the moment 
of inertia of the deckhouse and hull are A^, 1^, and 
respectively, and the distances of the respective centers of gravity 
from each other and from the deck are a, and o^a. 

At a distance z in the free body diagram of Figure 2, the moment 
and direct forces in the deckhouse and hull are M^-, , and M 0 , N 0 

respectively, with positive moments producing compression at the top of 
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the respective units. Direct forces and ^ are positive if they 
create tension. The external loads acting to the left of the section 
have a moment M. A shear force T of unknown magnitude will act on 



the underside of the deckhouse, and a similar force T will act in the 
opposite direction on the hull. Equilibrium of the portions of deck- 
house and hull in Figure 2, requires the relations: 



\ = -T, 


, M 1 = - a djT 


(5a) 


N 2 = T, 


M 2 = M - a cx 2 T 


(5b) 


Owing to the assumption that Navier's hypothesis is valid for the 


deckhouse and 


hull separately, the stresses can be determined at 


points 


at a distance 


or x^ from the respective center of gravities. 


In 


the deckhouse 3 


► 




a - Nl 
1_ A 1 


M i 


(6a) 


in the hull, 






N 2 

° 2 ' A 2 


M 2 

" I 2 X 2 


(6b) 



with tension stresses counted as positive. 

At the junction of deckhouse and hull, the longitudinal stress 
in the deck at this junction and the longitudinal stress in the deck- 
edge may differ; simple beam theory implies that they are the same. 
Shear lag in the deck plating may modify this using the following 
relation, 



( 3 ) 
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where r is the ratio of longitudinal deck stress at the junction to 
the stress at deck-edge. 

Furnishing Equations (5) and ( 6 ) with = -a = a 0 ^ and 

using Equation (3), the value of T is found to be: 



T = 



a a 2 Ii r 



1 l I 2 , I 1 I 2 . 2 , 2 , , _,2 T > 

T T l4,( “i rI i "iV 



M 



(7) 



T was defined as the total horizontal shear force acting between 
the left end of the deckhouse and the section at z, According to 
Equation (7), T is proportional to the moment M, and the unit horizontal 
shear (dT/dz) which will be transferred by rivets or welds from the hull 
to the deckhouse, will be 



dT 

dz 



a a 2 I 1 r 



I l I 2 , I 1 I 2 , 1,2 T , 2 T N 

“AT + "XT r + 3 (a 2 r h + a i I 2 ) 



V 



( 8 ) 



where V = 3 —, is the total shear force in the structure. 
dz 

After finding the value of T using Equation (7), and introducing 
this into Equations (5) and ( 6 ) , the longitudinal stresses O and O ^ 
at any point can be calculated easily. 
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CHAPTER II - GENERAL ANALYSIS OF TWO-CELL STRUCTURE 



Consider again the structure in Figure 1; differing from the 
treatment in the preceding section the assumption is made that any 
relative displacement of the deckhouse with respect to the hull will 
be resisted by the internal forces required to deflect bulkheads and 
transferse beams supporting the deckhouse. The deckhouse is considered 
as a beam elastically supported on the hull, and is further attached to 
the hull at deck level so as to enforce equal strains. External vertical 
loads and buoyancy will cause the structure to deflect, and this 
deflection can be described by the displacements y^ and y ^ of the center 
lines of the deckhouse and hull respectively, (Figure 3). In order to 
exclude motions of the entire vessel as a rigid body, y^ and y^ are 
defined as the relative displacements measured from a straight line 
C-D rigidly connected to the hull. As a result of this definition the 
displacement y ^ of the centroid of the hull at points C and D must 



always 


be zero. 

/ C.G. of Deckhouse 


r Deflected 
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~ 7 __ 
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FIGURE 3 
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It is assumed that the stiffness of bulkheads or deck beams, 

resisting relative vertical displacements of the deckhouse, is constant 

for the full length of the deckhouse, the magnitude of the stiffness 

being given by a spring constant K. K is defined by Bleich as being 

the force per unit length of deckhouse required to produce a relative 

deflection equal to one unit of length. Therefore the vertical 

# 

reaction between hull and deckhouse will be K(y^ - y^) per unit of 
length . 



The structure analysed here is shown in Figure 4. There are two 
beams having areas and A ^ and moments of inertia 1^ and and they 
are connected along C-D in such a way that both horizontal shear forces 
and vertical reactions can be transferred. It is assumed that, also, 
Navier’s hypothesis to be valid for the hull and for the deckhouse 
separately. 




FIGURE 4 



-24- 



This structure will be under the action of vertical loads 

on the deckhouse, on the hull (which includes buoyancy), shear forces 

Sqj and moments and M^. 

Using the n theorem of stationary potential energy 11 , the 

differential equations for the deflections y^ and y^ can be obtained. 

The total potential energy U consists of the internal strain energy V, 

♦ # 

and the potential energy U of the external forces. The total potential 
energy U = V + is"*" , 



1/2 






- 1/2 



[EI 1 yj2 + EI 2 y'^2 + El^a^ + r a^ 2 ) 2 + Uy 1 ~ yj 2 



+ E E* (y 1 ") 2 + E F*(y"') 2 + E G*(y"' y'") 
1 2 12 



- 2p 1 y 1 - 2p y ] dz + [M y']^ 2 - [S yj^ 2 
11 z z -H/2 * -i/2 



(9) 



U will be minimum if the variation 
6U = 0 



Using the rules of calculus of variations, the set of two simultaneous 
equations are derived. ^ 



See Appendix I, for the derivation of V and U^ 



2 



See Appendix II 
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VI 2 . IV EG* VI IV 

EE* y 1 1 + E^ + a x I A > y;L + K y ± - ~~ y^ + E r« 1 a 2 I A y 2 



- K y 2 = Pi 



(10. a) 



EG* VI , _ _ IV 

~T V 1 + E r a l a 2 h y l ' K y l " 



VI 2 2 , IV 

EF* y 2 + E(I 2 + r a 2 I A ) y 2 



+ K y 2 = P 2 



(10. b) 



The calculus of variations method also furnish the boundary 
conditions required to determine the arbitrary constants which will 
appear in the general solutions of the differential equations. 
Because of symmetry, there are six boundary conditions instead of 
twelve. For z = + 1/2 and z = - H/2 



?2 



= 0 



(11. a) 



2EE* y^ " + EG* y 2 ' ' = 0 



(11. b) 



EG* y| " + 2EF* y 2 " = 0 



» 1 I _ 



(11. c) 



TV 9 ft Itt II 

-EE* y^ v + E(I X + a* I A ) y x - y 2 + E r a i a 2 I A y 2 =0 (11. d) 

EG* IV 11 IV 2 . » i 

y/ + E r a x a 2 I A y x - EF* yj v + E(I 2 + r £ I A > y £ = -M 

(11. e) 



2 1 



V 2 » » > EG* V i i » 

EE* - E(I 1 a ± I A ) y x + — y 2 - E(r ^ « 2 I ft ) y 2 - 0 



<11. f) 
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II, A - NON-DIMENS IQNALIZATION OF THE EQUATIONS 

It is thought that employing non-dimensional equations and 
boundary conditions, the amount of algebra in the computation will 
be reduced and the results can be represented in a general form. 

Non-dimensionalization is made in the manner that is given 
by the following equations. 



y± = y ± l, i = 1,2 ' (12. a) 

z = z l (12. b) 

a = a H (12. c) 

A i = Jl 2 , i = 1,2 (12. d) 

I. = T7 £ A , i = 1,2 - (12. e) 

ii 

E£, i = 1,2 (12. f) 

M = M E l 3 (12. g) 

K = K E (12. h) 



The length of the deckhouse is selected as a characteristic length, 
because it is one of the most important parameters in the distribution 
of the longitudinal stresses. 

Substituting Equation (12) into Equations (10) and (11) , the 
following non-dimensional equations and boundary conditions are 
derived. 

-a i6 + a i4 + K y 1 - a y 2 + b y 2 - K y 2 = p x ( 13 * a > 

-a Yl + b ^ - Ky r a 26 ^ + a^ y 2 +.K ^ = ?2 (13.b) 
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Boundary conditions for z = + i/2 and z = - i/2, 



y 2 - o 

a l6 y l" 


+ a y’" = 0 




(14. a) 
(14. b) 


"fit , 
a y x + 


-?l» A 

a 26 y 2 ' 0 




(14. c) 


-IV 


. - 1 » -IV 


— IT 


(14. d) 


-a 16 y l 


+ a 14 y l - a y 2 


+ b y 2 =0 


-IV 


, -M “ -IV 


— ! ! - 


(14. e) 


-a y x + 


b y l ‘ a 26 y 2 + 


a 24 y 2 = - M 


-V 


-ITT -V 


- T ! T 


(14. f) 


a l6 y l ' 


a 14 y l + a y 2 ' 


b y ' 2 =o 



Expressions for non-dimensional coefficients for the equations and 
boundary conditions are given in the section entitled ’’Parameters 
Used for Computation”. 
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CHAPTER III - SOLUTION OF DIFFERENTIAL EQUATION FOR CONSTANT 

MOMENT M 



Considering the simple case that the loads p^, p^ and the 
shears and are zero, the only loads being M^ = M^ = M, 

Equations (13) are then homogeneous. 

Setting the determinant of the coefficients of the differential 
equations equal to zero, the charactericits equation will be derived, 
and it will be in the following form. 



4,8^ 6 ^ 4 2 v n 

r (r + a^ r + r + a^ r + a^) = 0 



(15) 



where, a^, a 2 > a^, and a^ are known constants. The roots of the 
characteristic equation is found for the models with different 
dimensions, and it is seen that the roots were always in the following 
manner. 

(16. a) 
(16. b) 
(16. c) 
(16. d) 



r l " r 2 ■ r 3 ’ r 4 " 0 
r 5 ‘ - y 3 
r 6 " ^3 



r 8 ’ \ 



r , - - - 1 r 2 



r 10 = - *1 + 1 ~<2 



r ll ■ >1 - 1 \ 



(16. e) 
(16. f) 

(16. g) 

(16. h) 



r l2 = Yi + 1 Y 2 



(16. i) 
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Then, keeping in mind that the problem considered is symmetrical 
with respect to the origin of the co-ordinate z, and using only 
symmetrical functions, the general symmetrical solution will contain 
only six arbitrary constants. This general solution is, 

y 1 = c 1 + c 2 z 2 + c^ sin y 2 z sinh z + c^ cos Y 2 z cosh Y^ 2 

+ c^ cosh 2 + c 5 cosh y^ 2 (17. a) 

y 2 = + c 2 z 2 + c^ sin Y 2 z sinh y^ + cos Y 2 z cosh y^ z 

- y^ c^ cosh Y 3 z - y^ c^ cosh y^ z (17. b) 



where. 



c 3 y i c 3 + y 2 c 4 



c 4 " y 2 C 3 + y i c 4 



and, y^, y 2 , y^, y^ are given in the following expressions, 



(a 16 h " a 14 f " K) 



- (a l6 m " a i4 g) 



-(am - bg) 

(-ah + bf - K) 



y i = 



(-ah + bf - K) 



(am - bg) 



-(am - bg) 

(-ah + bf - K) 



(18 .a) 
(18. b) 



(18. c) 



or, 
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(ah - bf + K) 
-(am - bg) 



( -a 26 h + a 2A f + K) 



" (a 26 m ‘ a 2A g) 



( " a 26 h + a 24 f + K) 



" (a 26 m ‘ a 24 g) 



(18. d) 



(a 26” ' a 24 g) 



<_a 26 h + a 24 f + K) 



(a l6 m ‘ a 14 g) 



<a 16 h - a 14 f - 



(-ah + bf - K) 



(am - bg) 



-(am - bg) 

(-ah + bf - K) 

-(am - bg) 

(-ah + bf - K) 



(18. e) 



or, 



(am - bg) 



• (a 26" - a 24 g) 



(ah - bf 4* K) 



(-a 26 h + a 2A f + K) 



(a 26 m " a 24 g) 



(_a 26 h + a 24 f + K) 



~ (a 26 m " a 24 g) 



( - a 26 h + a 24 f + K) 



(18. f) 



U 3 = 



(“ a l 6 Y 3 + a X 4 Y 3 + K) or (-a Y 3 + b Y 3 “ K) 



(-a yj + b yj - K) 



( - a 26 Y 3 + a 24 ^3 + « 



(18. g) 
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( ‘ a 16 Y 4 + a 14 Y 4 + K) or ( ' a Y 4 + b Y 4 - 

y = - = (18. 

(- a y 4 + b y 4 - K) (- Y 4 + Y 4 + K) 

Introduction of Equations (17) into the boundary conditions, 
leads to the equations to find the arbitrary constants given in 

t 

TABLE (1). 



TABLE 1 
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After finding the general coefficients, the stresses at any 
point along the deckhouse can be computed from the expressions for 
the moments M^, M^, and direct forces N^, , 



Mi = 


- E I x y” (19. a) 


m 2 - 


** E I 2 y 2 ’ # (19. b) 



Because, in this case p^ = p£ = 0. 



25 
h- 1 

11 


- N 2 * a ( “l’l tr “2 y 2 > (20> 


Here, it is 


important to notice that the solutions for and y 2 will 



be non-dimensional. 
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CHAPTER IV - SOLUTION OF DIFFERENTIAL EQUATION FOR 
FOR EQUALLY DISTRIBUTED LOADS 



In this section the case is considered of equally distributed 
loads and T> 2 > acting on deckhouse and hull, respectively, while 
the moments at the end of the deckhouse are = 0, Equilibrium 

requires external shear forces 

S C - - S D - f <»1 + V (21) 

at the ends C and D. The moment in the midship section due to the 
loads pi and p2 are 

^1 + ^2 2 

Mp = -^—2. I (22) 



M, 



a/2 


U* J 






^ Z 




1 
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ni i n 



% 



i i 1 1 b ” 

z. 



Equally distributed loads, 

Pi + p 2 



FIGURE 5 
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The loading being symmetrical, the general symmetrical solutions 
of Equation (13) are, 



-2 



y 1 = c^ + c 2 z + c 3 sin Y 2 z sinh y.^ z + c 4 cos Y 2 2 cosh y z 



+ c c cosh y_ z + c, cosh y. z 
5 3 b 4 



(Pi + P 2 ) -4 P-, 

z -f- 



24(a 14 + a 24 + 2b) 



(1 + v 5 )k 



(23. a) 



-2 



y 2 = c^ + c 2 z + c^ sin y 2 z sinh y^ z + c 4 cos y 2 z cosh y^ z 



- y 3 c^ cosh y 3 z - y 4 c^ cosh y 4 z 



(Pj + P2 } 



24(a l4 + a 24 + 2b) 



z 4 + 



5 P 2 



(i + p 5 )k 



(23. b) 



where c^, c^, y^, y£ » y^» and as given in the Equations (18) for 
the case of constant moment loading, and 






a 14 + b v 9 



5 a 24 + b v 10 



(24) 



But it should be mentioned that the value of "r" shear-lag factor will 
be different in this case, so all the coefficients will not, in 
general, have the same numerical value. 

Introducing general solutions Equations (23) into the boundary 
conditions Equations (14) will lead us to the same kind of equations 
given in Table 1. But, there will be some more terms, because of the 
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particular integral parts of the general solutions. These additional 
terms can be determined easily, using the particular integral parts 
of the general solutions. 

After finding the general coefficients, the further computation 
to find the stresses at any point along the deckhouse follows the 
pattern for the preceding section. The only difference being that 
the additional terms appear in the equations for and due to 
distributed loads p^ and p^ on the deckhouse and hull, respectively. 




( 25 ) 




' ( 26 ) 



( 27 ) 
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CHAPTER V - TOTAL SOLUTION OF DIFFERENTIAL EQUATION 

Considering the total loading for the system shown in Figure (5) , 
being equally distributed loads and acting on deckhouse and hull 

respectively, while the moments at the end of deckhouse are = M, 

it is possible to find the total solution 

The differential equations for this system are Equations (13) 

% 

and the boundary conditions are given by Equations (L4) for (+ — ) and 

i 

(- . The problem considered is symmetrical with respect to the 

origin of the co-ordinate z. The general symmetrical solutions of 
Equations (L3) can be given by Equations (23. a) and (23. b) for and 
Y 2 > respectively. The coefficients in the solutions (c^* c^, 

U^, y^) are given by Equations (18) and (24). 

But, it must be remembered that the values of the general 
coefficients (c-^, c^j c^» c^, c^, c^) will be different numerically 
than the values found for the two cases considered before. 

In the above approach to get the total solution, an additional 
assumption is made: Deck shear-lag factor "r 11 is considered to be 

constant over the length of the deckhouse and is given by Equation (4) , 
even though there is an applied constant bending moment at the end of 
deckhouse. 

Another approach to get the total solution is to make super- 
position to the solutions found for. constant moment and equally 
distributed loading cases. 

Total solution of differential equation is found by using both 
approachs explained above, and the results are compared in the 
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following chapter. 










D 
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CHAPTER VI - MODEL USED FOR COMPUTATION 



For this analysis a model was selected with the dimensions as 
shown in Figure (6). This model can be assumed as short deckhouse, so it 
will be possible to see more pronounced shear effects. Bulkheads are 
placed at equally spaced distances of 20 feet in the hull section. The 
thickness of the hull box girder plating is 0,5 inches and the thickness 
of the deckhouse plating and bulkheads is 0.25 inches. The material 
constants include a Young’s modulus of 30 x 10^, modulus of elasticity 
in shear of 11.5 x 10^ and a Poisson’s ratio of 0.3. 

The model is assumed to have a 15 foot draft with a corresponding 
hydrostatic upward distributed force of 17.143 tons/ft. The internal 
loading is arranged so as to provide equilibrium and a resultant 
symmetrical loading. Shear and moment diagrams for the total model are 
also provided in Figure (7). 

In a model with bulkheads, the main difficulty is the determina- 
tion of the deck stiffness or spring constant (K) . Since K was defined 
as the force per unit length of deckhouse required to produce a relative 
deflection equal to one unit of length, it is apparent that the value 
of K will, in reality, vary along the deckhouse length due to the 
emplacement of structural bulkheads. In order to simplify the use of 
the method, however, an average value of K must be determined. To 
achieve this end, the same approach* that was used in Reference (3) is 
followed. A symmetrical portion of the hull structure is modelled to 
include a bulkhead and attached deck and bottom plating, Figure (8). 
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^22.857 Tons/Ft 



D ^"11.429 Tons/Ft , 



:: i 



s_ 



/ 

17.143 Tons/Ft. 
5.714 Tons/Ft. 



I 



^ f_L 



5.714 Tons/Ft. 
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LOADING, SHEAR AND MOMENT DIAGRAM FOR 
MODEL WITH BULKHEADS 



FIGURE 7 
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SYMMETRICAL STRUCTURE FOR EVALUATION OF K 



FIGURE 8 
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This I-beam type structure is simply supported on its end and allowed 
to deflect under vertical forces applied to the hull-deckhouse con- 
nections. In the analytical procedure used for determining the forces 
needed to deflect the hull-deckhouse connections 1 inch, it must be 
kept in mind that although the deflection due to shear forces is 
negligible in most cases, in short deep metal beams the deflection 
caused by shear may become a significant portion of the total 
deflection. In this case shear contributes a major portion to the 
total deflection. 

4 

This analytic approach to K yielded a value of 1.58 x 10 psi 
for the present model under consideration. Sample calculations are 
provided in Appendix III. 

In Reference (3)>as a further check on K, a STRUDL program using 

*PSR f elements on the same model presented in Figure (8) was run. 

Arbitrary forces (F) were applied at the locations indicated, and the 

resulting deflection at the hull-deckhouse connection noted. The 

force was then scaled for a deflection of 1 inch. P was obtained by 

dividing F by the width of the flange (240 inches). Using Equation 

4 

K = 2P, K was found to have value of 1.53 x 10 psi. As it is said 
in Reference (3) , the disparity between the STRUDL K and the analytical 
K can be attributed to the fact that in the analytical approach the 
deflection calculations apply to the neutral axis of the beam only. 

The value of K that is used in the computations was the value 
found by employing STRUDL program. 
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VI, A - PRESENTATION OF THE RESULTS 

Presentation of the stress distributions for Bleich f s method 
and for the method developed are shown in Figures (9), (10), (11), (12), 
(13), (14), and (15) on the following pages. Figures (9) and (10) 
show the comparison of longitudinal stress distribution for constant 
moment loading at midship and at the end of deckhouse, respectively. 
The comparison of the results for equally distributed loading are shown 
in Figures (11) and (12) at midship and at the end of the deckhouse, 
respectively. Figures (13), (14), and (15) show the comparison of 
stress distributions for the total solutions. 

For the computation of the results, the model shown in Figure (6) 
is used for both methods. 

"ACCESS II" Primer Operations in Linear Algebra for the Interdata 
Computer in Joint Computer Facility is used to find the roots of the 
characteristic equations and to solve the general coefficient matrices. 

The explanation about Bleich's method can be found in References 



(1) and (3). 
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"CONSTANT MOMENT LOADING" 

COMPARISON OF LONGITUDINAL STRESS DISTRIBUTION AT MIDSHIP 




COMPRESSION TENSION 



All Stresses are in psi x 10 



2 



FIGURE 9 



RESULTS 

BLEICH'S 

RESULTS 



-46- 

"CONSTANT MOMENT LOADING" 



COMPARISON OF LONGITUDINAL STRESS DISTRIBUTION AT THE END OF DECKHOUSE 
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All stresses are in psi x 10 

FIGURE 10 
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"EQUALLY DISTRIBUTED LOADING" 

COMPARISON OF LONGITUDINAL STRESS DISTRIBUTION AT MIDSHIP 




All stresses are in psi x 10 
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RESULTS 



FIGURE 11 
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"EQUALLY DISTRIBUTED LOADING" 

COMPARISON OF LONGITUDINAL STRESS DISTRIBUTION AT THE END OF DECKHOUSE 




COMPRESSION TENSION 



All Stresses are in psi x 10" 



FIGURE 12 
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" TOTAL SOLUTION" 

COMPARISON OF LONGITUDINAL STRESS DISTRIBUTION AT MIDSHIP 




FIGURE 13 
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"TOTAL SOLUTION" 



COMPARISON OF THE RESULTS AT MIDSHIP 




COMPRESSION TENSION 



All stresses are in psi x 10 
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FIGURE 14 
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"TOTAL SOLUTION" 

COMPARISON OF LONGITUDINAL STRESS DISTRIBUTION AT THE END OF DECKHOUSE 
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All stresses are in psi x 10 
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FIGURE 15 
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VI. B - COMPARISON OF THE METHODS 



In Figure (9), the longitudinal stress distributions at midship 
found by using two methods were plotted. In this figure, it can be 
easily seen that both results are almost the same except 400 psi 
difference in stresses at deckhouse top. The difference is 50 psi at 

9 

main deck level. 

Figure (10) shows the comparison of two methods at the end of 
deckhouse. In this case, the difference in stresses at deck level is 
320 psi. At deckhouse top, there is a big gap between the values of 
stresses found by using two methods. At hull bottom, the difference 
is about 30 psi. But, the value of stress found by using the theory 
developed is smaller than the value found by Bleich f s method. 

Figures (11) and (12) show the comparison of longitudinal 
stresses for equally distributed load. Looking at the results shown 
by solid line in Figure (11) , shear-lag effect at main deck level can 
be seen easily. All the values for stresses are smaller than the 
results found by Bleich's method. But, the differences are not so big 
at mid ship. 

In Figure (12), it is not difficult to recognize the big 
differences between the results found by using both methods. The 
differences are approximately 1290 psi, at deckhouse top, and 690 
psi at hull bottom. At main deck level, there is about 310 psi and 
400 psi difference between the results at hullside and at the junction 
of main deck and deckhouse, respectively. 
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Figure (13) shows the comparison of the total solution at 
midship. The strudl results were taken from Reference (3). The 
agreement among the results is reasonable. The stress distributions 
were transformed into moments and checked against equilibrium condition; 
that is, the values of the moments obtained corresponded to the value 
on the bending moment diagram, except for strudl results and for the 
total solution using r = 0.76, strudl results give 9% bigger than the 
moment at midship, and the total solution with r = 0.76 gives 6% smaller 
bending moment. This was expected, because shear-lag factor r = 0.76 
is not true for the whole system with applied constant bending moment 
at the end of deckhouse. 

In Figure (14), the total solutions found by using r = 0.76, 
r - 1.0 and the total solution by superposition were compared. In 
this figure, it is seen that the total solution found by using r = 1.0 
takes the average of the stresses at deck level. 

Comparison of the total solutions at the end of deckhouse were 
shown in Figure (15) . The values of the stresses found by using the 
theory developed were always smaller than Bleich’s results. 

As a conclusion of the comparison of the methods, inclusion of 
the shear effects into two beam theory did not change the values of 
the stresses at midship. But, at the end of deckhouse, more 
distinguishable differences were found. 
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CONCLUSIONS AND RECOMMENDATIONS 



Even though the method which includes shear effects requires 
more elaborate work, it is seen that it is possible to use it for 
computing the longitudinal stresses at any point along the deckhouse. 
There seems to be no reason to ignore these effects in design procedures 
for short and moderate deckhouses where the shear effects could be 
pronounced . 

After examining the comparison of the results found by using 
both methods at midship, it is possible to conclude that there is some 
indication to use Bleich T s method for design purposes in the computa- 
tion of longitudinal stresses at midship. The simplicity in his compu- 
tation is a prime factor. But, for the stress solution at the end of 
the deckhouse where more distinguishable differences were found by 
employing both methods, the inclusion of shear effects into two-beam 
theory being used is necessary. The results may become more realistic. 

In the theory developed, the shear-lag effects in the deck 
represented by "r" is important only in the case of equally distributed 
load. Application of Equation (4) is a very simplified means of 
estimating "r", because it was determined by a box-girder analysis of 
the hull alone in Reference (2) . Even though it may require a much 
more elaborate analysis, it is recommended to improve on the estimation 
of the shear-lag parameter. 

In this analysis, the stiffness of bulkheads or deck beams 
resisting relative vertical displacements of the deckhouse was assumed 
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constant for the full-length of the deckhouse. It is known that the 
deck stiffness K vary along the deckhouse length due to the presence 
of structural transverse bulkheads. The inclusion of this variation 
in K in the two-beam deckhouse theory needs further investigation. 

Finally, it is recommended to built a physical model similar to 
the one considered in this theoretical analysis. Then the strain and 
deflection measurement under similar loading conditions could be 
employed for verification of the results. 
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APPENDIX - I 

I. A - DERIVATION OF STRAIN ENERGY OF STRUCTURE 

Denoting by £^ and £ ^ average longitudinal strain in the 
deckhouse and hull, respectively, the strain energy of the longitudinal 
stresses will be as follows: 

Deckhouse: 



there will be energy stored in the bulkheads or deck beams which resist 
the relative vertical displacements of deckhouse and hull; this part 
of the strain energy can be expressed by the spring constant K in the 
form 



The strain energy due to the shear deflection caused by vertical 



%/2 



» 



E 

2 



(A i e i + h y i' 2) dz 



M ,/ 2 



Hull: 




E 

2 



( A 2 <=2 + X 2 y ’ ' 2) dz 



2 



1/2 



The strains are counted positive if they represent elongation. 



In addition to the strain energy of the longitudinal stresses, 



1/2 



1 

2 



K( y i - y 2 ) 2 dz 



U/2 



shear forces will be as follows: 



Deckhouse: 



1 

2 
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Hull: 



1 

2 



i/2 . * 2 

k 2 V 2 
a 2 g 



dz 



-Jl/2 



where, k* and k* are shear deflection constants and for wide flange 
box girders, can be assumed that they are equal to 1.0. 

The total strain energy V is 



V - I 

V 2 



i/2 



t \ 2 t i |2 , . 2_ 1 * 2 , K / 

lA l S 1 + *1 V + A 2 C 2 * h H + I < y l - y 2 ) 



-i/2 



1 2 , 1 2 . . 

+ 2 EG A X V 2 + 2EG A £ V 2 ] dz 



(a) 



The stresses in the deckhouse and hull can be expressed by 
the average strains and anc * by the second derivatives of the 
deflections y| ? and y^ ? 

Deckhouse: 0 = Ee^ + E y| ? x^ 

Hull: 02 - E02 + E y” X 2 

Using the relationship given in Equation (3) in Chapter I 
for the longitudinal stress in the deck at the junction, and the 
longitudinal stress in the deck-edge, 

£ 1 - r e 2 = aa i *1 + r 3 a 2- y 2 f (b) 

Further, the longitudinal resultant of all stresses in the deckhouse 
N^, must be equal to EA^ e^, and, similarly, n 2 = ea 2 e 



The 
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resultant of all longitudinal forces, + N 2> must vanish. 



N x + N 2 = E(A X + A 2 e 2 ) *= 0 



(c) 



By means of Equations (a) and (b) , and e 2 can be expressed as 
follows: 



: 1 ' rA, 2 + A, (C, 1 y i' +r °2 y 2 ') 



(d) 



a A 



e 2 ” rA, + A„ ^1 J 1 



1 , » f » f X 

(a n y, + r a 2 y 2 ) 



(e) 



1 2 

Substituting expressions for E^ and e 2> and and v 2 (Equations 
1, and m, in I-C) into the Equation (a) for total strain energy V, 

f Jt/2 



V = E 
2 



[I 1 y l 2 + T 2 y 2 2 + V a i y l' + r “2 y 2 ,)2 + f (y l " y 2 )2 



-a/2 



. * » 1 2 . » » » ^ TfT -. 

+ E* y^ + F* y 2 + G* y 2 ] dz 



f » f 2 



it? tit 



(f) 
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I.B - POTENTIAL ENERGY U OF EXTERNAL FORCES 

w 



Counting and p 2 positive if acting downward, their potential 
energy is, 

rl/2 

( p i + p 2 y 2 ) dz 

-in 

The shear forces S and and the moments M and >L act 

c D cD 

immediately outside points C and D, and their potential energy will 
depend on the vertical displacements y^ c and y^, and on the rotations 
of* the end surfaces of the hull, y^, an< ^ ^2D* 

Taking into account the direction of. the shears and moments 
shown in Figure 6 in Chapter II, the potential energy will be; 



' S D y 2D + S c y 2c + “d y 



. 1/2 



., 1/2 



2D M c y 2c " ^ Sy 2^1/2 + 2^1/2 



and the potential energy of the external load is; 



, n 1/2 



. 1/2 



fl/2 



»« ■ ^2h,2 - lS ?2h/2 - (pl y l + p2 y 2 ) dz 

- 1/2 



(g) 
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I,C - EXPRESSION FOR LONGITUDINAL FORCES AND 
HORIZONTAL AND VERTICAL SHEARS 



The resultants N^ and N^ of the longitudinal stresses in the 
deckhouse and hull are N^ = EA^ and N 2 = EA^ £2* an< ^ using 
Equations (d) and (e) 



ET 

w i ■ — <o i * 1 ' + r “2 ^'> 



ei a 

N 2 = " T" (CX 1 y l' + r a 2 y 2 ' } 



where, I - a 



2 A 1 A 2 



(r A x + A 2 ) 



(h) 

(i) 



T being the total shear force from the left end of the deckhouse to a 
point having the co-ordinate z, equilibrium requires, 



El 



T = - N x = 



A / t 1 

- («! y x 



+ r a, 



» » V 

) 



2 J 2 



and the unit horizontal shear will be, 



dT 

dz 



El 



A / » » 1 , » » * v 

- ( 0l y x + r a 2 y 2 ) 



(j) 



(k) 
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FIGURE A 



To obtain the expression for the vertical shear V^, in the 
deckhouse, considering an element of the deckhouse of length dz , as in 
above Figure A. 

Equilibrium of moments with respect to the centroid requires, 



V 



1 



a a 



dT 
1 dz 



dMj^ 

dz 



Substituting — from the above Equation (k) and using 



M 1 " EI 1 (y l + A G 

lw 



V 1 - - E(I 1 + a l I A )y” ' -Oj^r EI a yj" 



( 1 ) 
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The vertical shear force in the hull, considering Figure B below, 
similarly: 

V 2 = - a L a 2 EI a yj” - E(I 2 + r a* l A ) y'" (m) 




M„ 




dz 



FIGURE B 

In the derivation of the Equations (1) and (m) , attention must 

dM. 

be paid to the following point; when the derivative of — — (for i = 

1,2) is taken, the contribution from the second term in the moment 

equation will be zero in both cases. Because, in the following 

chapters, the cases considered for solving the equations are constant 

moment loading (p = p 0 = 0), and equally distributed loading 

dp. dp 1 2 

( dT- IT' 0) - 
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APPENDIX II - APPLICATION OF CALCULUS OF VARIATIONS 



The problem considered was in the following manner, 



U = 6 



1/2 



F(z, y x , yj, yj\ yj", y 2> y' 2 , y' 2 ' , y' 2 ") dz 



- 1/2 



+ w y'llm - I s **1 m 



To get the set of two simultaneous differential equations, 



Euler-lagrange equation is employed for and in the following 



manner , 



Fy - dx + ~~2 *y' -73 = 0 

dx dx 



where , 



Fy = 



3F 

V 



3F 

p y = gT. et c. 



To get the natural boundary condition equations, the following 
equations are used for y^ and y 2 



U" = 0 

. r 3F d , 3F A . 

n 3y' ' dx ^3y ' ' ' " 0 

„ r 3F d , 3F N , d 2 , 3F N1 
n[ 3y* " dx ( 3y ' ' + dx 2 ( 3y ,,,)] “ 0 



APPENDIX III - CALCULATIONS FOR DETERMINING K 



The following is an analytical appraoch to the determination 
of the spring constant (K) for a model with bulkheads spaced 20 feet 
apart. 

Basic Nomenclature 
A Area of beam cross section 

K 1 Factor depending on shape of beam cross section 

p Distributed load at hull-deckhouse connection 

V Vertical shear due to actual forces 

v Vertical shear due to load of one pound acting at the section 

where the deflection is to be determined 
Y_, Deflection due to internal moments 

Y Deflection due to shear 
T 

Y t Total deflection 

x Distance along length of beam 

< > Indicates singularity functions 

The equation for the deflection due to the internal moments 
(Y ) is calculated through the use of singularity functions. The 
expression for the deflection due to shear (Y^) is obtained through 
the use of the method of unit loads as described in Reference (11) . 

The total deflection (Y T ) is then expressed as the sum of Y^ and Y 
and set equal to 1 inch. The equation is solved for p (the distributed 
load acting on the hull-deckhouse connections) . The expression for 
K is twice the value of p. 
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Deflection Due to Moment Only 



d 2 Y , i 

El — = 240 [p - p <x - 120> - p <x - 360> ] 

dx X 



EIY M " 240 Ir - ^ - "6^ P>3 ' + c, x + c 2 ] 



B.C Y„ = 0 when x = 0, 480" 

E = 30 x 10^ psi 
1 = 8,748,005 in 4 

Y = - (2.105 x 10 ^)p at x = 120" 
M 

Deflection Due to Shear 



Y 

T 



1 

K’ 



V 

— dx (method of unit loads) 



v = I p(240)x 
T K' AG 



K’ = 



10(1 + y) (1 + 3m) 

6(2 + 12m + 25m 2 + 15m 3 ) + (11 + 66m + 135m 2 + 90m 3 ) + 



3 Oran 



i 2 (l + m) +5 + 9m) 



where 




I-BEAM CROSS SECTION 
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t- = 0.5"; t = 0.25"; h = 360"; b = 240" 
f w 

m = 2.66; n = 0.666 
K f = 0.262 

It is customary to assume that only webs of structural shapes, such as 
channels and I beams resist shearing stresses* because shear stresses 
are small in flanges. 

Y - - (10.55 x 10" 5 )p 

. y t - y m + y t - lH 

1" = - (12.655 X 10 _5 )p 
p = 0.7902 x 10 4 
K = 2p 

K = (1.58 x 10 4 )psi 
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APPENDIX IV - SAMPLE CALCULATIONS 
IV. A - FOR CONSTANT MOMENT LOADING 
a = 306", A x = 150 in 2 , A 2 = 840 in 2 

I 1 = 534,600 in 4 , I 2 = 19,440,000 in 4 , I = 11,917,309 in 4 
I = 31,891,909 in 4 

7 ' 

M = M = 19.2 x 10 in-lbs 
P 

K = 1.53 x 10 4 

= 0.4118, a 2 = 0.5882, r = 1.0 (shear-lag factor) 



b = (1.392) 


x 10 6 , 


a l4 = (1.232) 


x 10 6 , 


a 2A = (1.136) x 10 5 


a = (5.685) 


x 10 -8 , 


a l6 = (2.335) 


x 10" 8 , 


a 26 = (3.129) x 10 -7 



4,8 ± 6 4 2 ^ n 

r (r + a^ r + a 2 r + a^ r + a^) = 0 

a^ = -(1.208) x 10 2 , a 2 = (2.959) x 10 3 , a^ = -(5.630) x 10 4 
a Q = (1.924) x 10 6 

Roots of the characteristic equation are as follows: 

r l = r 2 = r 3 = r 4 = ° 

r 5 = " r 6 = " Y 3 = ' 5,843 

r 7 = - r 8 = - y A = - 9.656 

r 9 = - Yl - i Y 2 = - 3.266 - i 3.731 

r 1Q = - y 1 + i y 2 = “ 3.266 + i 3.731 

r u - Y x - i Y 2 = 3.266 - i 3.731 

r 12 = Y 1 + 1 Y 2 = 3,266 + 1 3,731 
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f = -(583.35), g = (158.58), h = 5762.89, m = 13,700.88 
n 1 = (-) 3.253, n 2 = 24.370, r ± = (-) 101.55, r 2 = (-) 67.45 
= 2236.18, S 2 = 2694.43 

y 1 = (-) 0.20974, y 2 = (-) 0.00449, y 3 = (-) 0.88467 
y 4 = 0.22336 

t x = (-) 1.166 x 10" 6 , t 2 - (-) 7.449 x 10" 7 , t 3 = 7.982 x 10 

t. = 7.359 x 10 -7 , t. = 3.413 x 10" 6 , t, = (-) 39.09 x 10 _6 

H 3 O 

t ? = (-) 3.371 x 10" 6 , tg = 24.899 x 10~ 6 , t g = 12.216 x 10 _5 
= 9 . 302 x 10 5 

v 1 = 7.364 x 10” 8 , v 2 = 1.065 x 10” 8 , v 3 = 33.375 x 10" 8 

v. = (-)(1.305) x 10 8 , v = 2.463 x 10 8 , v £ = 0.921 x 10 8 

= 11.444 x 10 8 , Vg = (-) 1.146 x 10 8 , = 2.624 x 10 8 

v^q = 12.755 x 10 8 

For the B.C's equations; (z = 0.5) 

sin sinh y^z = 2.34891 
cos y 2 z cosh y^z = (-) 0.76987 
sin y^z cosh y^z = 2.53634 
cos Y 2 z sinh y^z = (-) 0.71298 
sinh Y 3 z = 9.257 
cosh Y 3 Z = 9.404 
sinh y^z = 62.540 
cosh Y/ + z = 62.548 
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Substituting everything into the B C's equations shown in Table 1, 
the following matrix is derived, 



— 




















1.0 


0.250 


-0.49611 


0.1509 


8.3199 


-13.970 




C 1 


= 


0.0 


0.0 


0.0 


-3.5181 


-1.049 


136.010 


599.812 




C 2 




0.0 


0.0 


0.0 


2.549 


1.297 


616 . 40 


-735.20 




C 3 




0.0 


0.0 


5.249 


-22.08 


-94.45 


-16.22 


-418.6 




* 

C 4 




0.0 


0.0 


25.51 


11.25 


56.73 


16.12 


416.3 




C 5 




-3.704 x 10 -i 


0.0 


0.0 


37.62 


14.88 


9.321 


769.2 




C 6 




0.0 























Results from the solution of above matrix, 

C;L = 3.494 x 10" 5 

c 2 = - 1.290 x lO -4 

c 3 = 3.070 x lO -6 

c = - 7.899 x 10 -6 
4 

c 5 = 7.021 x lo" 9 

c, = 2.549 x 10 -9 
o 



The value of is not used in the calculations, as it was pointed 
out by Bleich in Reference (1) ; it describes only a rigid-body motion 
of the structure not required for the purpose of this analysis. 
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For z = 0 (Midship) 



•r^ = - 1481.439 psi 
A 1 

M 1 

— x 1 = (-) 494.586 psi; x.. = - 126.01" (main-deck level) 

h 

M 1 

— x^ = + 211.948 psi; x^ = 54.0" (deckhouse top). 



. „ 11 

so; using a. = — — x 

A 1 1 

= - 1481.439 + 494.586 = - 986.85 psi (main-deck level) 
0 1 = - 1481.439 - 211.948 = - 1693.38 psi (deckhouse top) 



= 264.542 psi 

A 2 

m 2 

— x_ = 1251.396 psi; x„ = 179.989 (main-deck level) 
L 2 Z 

^2 ii 

— x 2 = _ 1251.396 psi; x 2 = ~ -*-79.989 (hull-bottom) 



Using, o 2 





x 



2 



0 2 = 264.542 - 1251.396 = - 986.54 psi (main-deck level) 
0 2 = 264.542 + 1251.396 = 1515.938 psi (hull-bottom) 
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For z = 0.5 (At the end of deckhouse) 



■— « - 982.177 psi 
A 1 



M.. ii 

— x^ = 348.770 psi; = - 126.01 (main-deck level) 



M 

— x^ = - 149.461 psi; x^ = 54" (deckhouse top) 



= - 982.177 - 348.770 = - 1330.94 psi (main deck level) 
a 1 = - 982.177 + 149.461 = - 832.71 psi (deckhouse top) 



—■ = 175.388 psi 
A 2 

M 

— x 2 = 1506.336 psi; x^ = 179.989 (main-deck level) 



2 " 

— x 2 = - 1506.336 psi; x 2 = “ 179.989 (hull-bottom) 

O 2 = 175.388 - 1506.336 = - 1330.948 psi (main deck level) 



O' 2 - 175.388 + 1506.336 = 1681.724 psi (hull-bottom) 
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IV. B - FOR EQUALLY DISTRIBUTED LOADING 

All the values of parameter for mathematical model will be 
the same except, 

r = 0.76 (shear-lag factor) 

2 A 1 A 2 7 A 

1 , = a V ■ = 1.23670 x 10 in 

A rA^ i ^2 * 

I = I ± + I 2 + I = 32,341,600 in 4 

b = 1.09790 x 10~ 6 , a l4 - 1.26918 x 10 -6 , a 24 = 1.09432 x 10 

a = 5.28018 x 10" 8 , & = 2.48364 x 10" 8 , a 26 = 2.82865 x 10 

Coefficients in the characteristic equation, 

a 3 = - 1.21261 x 10 2 , a 2 = 2.98730 x 10 3 

a 1 = - 4.96456 x 10 4 , a Q = 1.730691 x 10 6 

Roots of the characteristic equation are as follows. 

r l = r 2 = r 3 = r 4 - ° 

r 5 = - r 6 = - = - 5.882 

r 7 = " r 8 = “ = - 9.624 

r 9 = - Y x - i Y 2 = " 3,183 " 1 3 - 621 

r^Q = - Y-^ + 1 Y 2 = “ 3.183 + i 3.621 

r u = Y x " 1 Y 2 = 3 - 183 " 1 3.621 

r 12 = Yj + i Y 2 = 3,183 + 1 3 -621 
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f = - 522.480, g = 137.393, h = 4724.149 
m = 11,634.397 



3 

H 1 

II 


- 2.980, n 2 = 23.051 


r i = 


- 92.955, r 2 = - 62.581 


cn 
1— 1 

II 


1895.044, S 2 = 2329.223 


II 

1 — 1 


- 0.20766, V 2 = ~ 0.01375, = - 0.72377 


11 


0.25226 


fc l = 


- 1.289 x 10' 6 , t 2 = - 8.006 x 10 -7 , t 3 = 3.085 x 10~ 7 


S = 


7.331 x 10 -7 , t c = 3.697 x 10 -6 , t, » - 26.240 x 10 _7 

D 0 


t 7 = 


- 3.617 x 10" 6 , t 0 = 26.307 x 10" 6 , t_ = 12.570 x 10 _5 

0 y 


rr 
1— 1 
O 


= 10.272 x 10 -5 


v i = 


6.305 x 10~ 8 , v 2 = 1.151 x 10 -8 , v 3 = 25.753 x lO -8 


C 

-O 

II 


-1.855 x 10 -8 v 5 = 2.063 x 10~ 6 , v & = 0.991 x 10~ 6 


V 7 ~ 


- 9.018 x 10" 6 , v 8 = - 1.667 x 10~ 6 , v g = 2.367 x lO -6 



v 1Q = 12.041 x 10 

For the B.C's equations; (z = 0.5) 
sin y 2 2 sinh z = 2.2854 

cos y 2 z cosh Y 1 z s • 0.6060 

sin Y 2 z cosh y^ z = 2.4832 
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cosh Y 2 z sinh Y-l z = “ 0.5578 
sinh Y 3 z = 9.43150 
cosh Y 3 z = 9.48439 
sinh Y^ z = 61.517 
cosh Y 4 z = 61.526 



Substituting everything into the B C’s equations shown in 
Table 1, the following matrix is derived. 



1.0 0.250 -0.4829 0.0944 



6.864 -15.55 



0.0 


0.0 


-3.648 


-1.268 


121.0 


631.2 


0.0 


0.0 


11.75 


16.48 


4943.0 


-10170. 


0.0 


4.734 


-7.452 


-62.61 


-38.86 


-427.8 


0.0 


24.08 


7.676 


62.31 


35.38 


291.5 


0.0 


0.0 


36.94 


18.50 


22.73 


411.7 



c 



/ 



Results from the solution of above matrix, 

,-5 



= 5.678 x 10 
C 2 = - 1.000 x 10 



-4 



c^ = 7 . 168 x 10 



-6 



c 4 = -1.015 x 10 



-6 



Cj. = -6.462 x 10 



-7 



c, = 3.548 x 10 
o 



-8 



2.322 x 10 



-5 



-8.066 x 
-3.487 x 

-4.539 x 
-2.430 x 
2.459 x 



10 

10 



-5 

-3 



10 



-4 



10 



-3 



10 



-4 
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The value of is not included in the computations. As explained 
before, it describes only a rigid-body motion of the structure, not 
required for the purposes of this analysis. 

For z = 0 (Midship) 

N 1 

~ ■ - 957.677 psi 
A 1 

M. 

Y~ = - 200.573 psi; = - 126.01 (main deck level) 



M 1 

— x- = 85.592 psi; x. = 54" (deckhouse top) 

1 

using, o x - — - — 

a 1 = - 957.677 + 200.573 = - 759.103 psi (main deck level) 
0^ = - 957.677 - 85.592 = - 1045.62 psi (deckhouse top) 

N 2 

— = + 171.371 psi 
A 2 

M 

— = 1224,315 psi; = 179,989 (main deck level) 



2 " 

~ x 2 = - 1224,315 psi; x^ = - 179.989 (hull-bottom) 




O = 171,371 - 1224.315 = - 1052.94 psi (main deck level) 
O = 171.371 + 1224.315 = 1395.686 psi (hull-bottom) 
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For z = 0.5 (At the end of deckhouse) 

N 1 

— = - 520.778 psi 
A 1 

M 

— = - 148.120 psi; x^ = - 126.01 (main deck level) 

M 

— x. = 63.474 psi; x.. = 54" (deckhouse top) 

1 

= -520.778 + 148.120 = - 372.658 psi (main deck level) 
= - 520.778 - 63.474 = - 584.252 psi (deckhouse top) 

!t 

179.989 (main deck level) 

II 

= - 179.989 (hull-bottom) 

o 2 = 92.996 - 583.474 = - 490.478 psi (main-deck level) 

©2 = 92.996 + 583.474 = 676.470 psi (hull-bottom) 



N„ 



M„ 



= 92.996 psi 



~ x 2 = 583.474 psi; x 2 = 



M 2 

— x 2 = ~ 583.474 psi; x^ 

n h h 

°2 ' A 2 I X 2 
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IV. C - FOR TOTAL LOADING 



All the values of parameter for mathematical model will be the 
same except, 

r = 0.76 (shear-lag factor) 

I. = 1.23670 x 10^ in 4 

A . 

I = l 1 + I 2 + I = 32,341,600 in 4 

The other coefficients and parameters will be the same as for 
equally distributed loading case. 

Substituting everything into the B.C's equations, the following 
matrix will be derived 



S' 

1.0 


0.250 - 0.4829 


0.0944 


6.864 


- 15.15 




S 

C 1 




/ 

2.322 


X 


10 -*' 


0.0 


0.0 - 3.648 


- 1.268 


121.0 


631.2 




c 2 




- 8.066 


X 


10 -5 


0.0 


0.0 11.75 


16.48 


4943.0 


- 10170.0 




C 3 




- 3.487 


X 


io " 3 


0.0 


4.734 - 7.452 • 


- 62.21 


- 38.86 


- 427.8 




C 4 




- 4.539 


X 


io " 4 


0.0 


24.08 7.676 


62.31 


35.38 


291.5 




C 5 




- 6.134 


X 


io -3 


0.0 


0.0 36.94 


18.50 


22.73 


411.7 

> 




c 6 




2.459 

k 


X 


io -4 



Results from the solution of above matrix, 
c., = 9.244 x 10 _5 



c 2 = ~ 2.286 x 10 
Cg = 1.234 x 10 3 
c. - - 1.145 x 10 
Cj. = - 6.169 x 10 
c 6 = 3.878 x 10 -8 
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For 2=0 (Midship) 

h 

= - 2049.499 psi 
M 

— x, = - 492.638 psi; x. = - 126.01 (main deck level) 

I 1 

M 1 • . 

— - x. = + 211.113 psi; x.. = 54" (deckhouse top) 

I 1 

0 1 = - 2049.499 + 492.638 = - 1556.861 psi (main deck level) 

0 ^ = - 2049.499 - 211.113 = - 2260.61 psi (deckhouse top) 

N 2 

~ = 365.982 psi 
A 2 

m 2 

— x„ = 2414.481 psi; x 2 = + 179.989 (main deck level) 
l 2 

M 2 

— x 2 = - 3414.481 psi; x 2 ~ ~ 1^9.989 (hull-bottom) 

0 2 = 365.982 - 2414.481 = - 2048.49 psi (main deck level) 

o 2 = 365.982 + 2414.481 = 2780.46 psi (hull-bottom) 

For z = 0.5 (At the end of deckhouse) 



(main deck level) 
(deckhouse top) 



—■ « - 1035.029 psi 
A 1 



11 

-p x x = 463.73 psi; = 126.01 



x = - 198.725 psi; x = 54" 
1 1 
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0 1 = - 1035.029 - 463.730 = - 1498.759 psi (main deck level) 
0 = - 1035.029 + 198.725 = - 836.30 psi (deckhouse top) 



— = 184.826 psi 
A 2 



M 

— x 2 = 2156.883 psi; x ^ = 179.989 (main deck level) 

M 

— x^ = - 2156.883 psi; = 179.989 (hull-bottom) 

O ^ = 184.826 - 2156.883 = - 1972.05 psi (main deck level) 



= 184.826 + 2156.883 = 2341.709 psi (hull-bottom) 
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